Harmonic Generation from Laser-Irradiated Clusters 



o 
o 

u 

Or 
< 

in 

(N 

5 1 

a* 



. 

11: 

o ■ 

& : 

Ph. 



O 

o 
o 



X 



M. Kundu, 1 S.V. Popruzhenko, 1,2 and D. Bauer 1 

1 Max- Planck- Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 
^Moscow State Engineering Physics Institute, Kashirskoe Shosse 31, 115409, Moscow, Russia 

(Dated: February 5, 2008) 

The harmonic emission from cluster nanoplasmas subject to short, intense infrared laser pulses 
is analyzed by means of particle-in-cell simulations. A pronounced resonant enhancement of the 
low-order harmonic yields is found when the Mie plasma frequency of the ionizing and expanding 
cluster resonates with the respective harmonic frequency. We show that a strong, nonlinear resonant 
coupling of the cluster electrons with the laser field inhibits coherent electron motion, suppressing 
the emitted radiation and restricting the spectrum to only low-order harmonics. A pump-probe 
scheme is suggested to monitor the ionization dynamics of the expanding clusters. 
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During the last decade, the study of rare gas and metal 
clusters interacting with intense infrared, optical and ul- 
traviolet laser pulses has emerged as a promising research 
field (see [l| for a state of the art-review). Genera- 
tion of fast electrons and ions, production of high charge 
states, generation of x-rays and fusion in deuterium clus- 
ters have been observed in laser-cluster experiments with 
laser intensities up to 10 20 W/cm 2 . The hot and dense 
nonstationary, nanometer-sized plasmas far from equi- 
librium whose dynamics occur on the femtosecond time 
scale — so-called nanoplasmas — are new physical objects 
with interesting properties very different from plasmas 
where finite-size effects can be neglected. The probably 
most important feature of these nanoplasmas is the un- 
surpassed efficient energy transfer from the laser light to 
the charged particles [2| . 

While the energy absorption by cluster nanoplasmas 
has been widely studied both in experiments and theory, 
much less attention has been paid to the laser harmonic 
emission from such systems. Meanwhile, it is known 
from the physics of intense laser- atom interaction that 
the effects of multiphoton absorption, leading to so-called 
above threshold ionization and high order harmonic gen- 
eration, are intimately related and can be described as 
different channels of the highly nonlinear laser-atom cou- 
pling. In intense fields, the laser-cluster coupling is also 
known to be highly nonlinear. In fact, nonlinear reso- 
nance [H has been shown (both in simulations and simple 
analytical models 0, HI, 0, 0]) to be of particular impor- 
tance for the energy transfer from the laser pulse to the 
electrons of the nanoplasma and the subsequent outer 
ionization. Most naturally the question arises whether 
laser-driven clusters can be an efficient source of high- 
order harmonics as well. Up to now, very few experi- 
mental results on harmonic generation (HG ) fr om cluster 
targets have been published. In Refs. HG from 

rare-gas clusters irradiated by infrared pulses of moder- 
ate intensity (~ 10 13 — 10 14 W/cm 2 ) was measured. It 
was shown that under such conditions harmonics can be 
generated up to higher orders and with higher saturation 



intensities than in a gas jet. However, the applied inten- 
sities were definitely not high enough to create a dense 
nanoplasma inside clusters so that the effects observed in 
Refs. H 0, U3] should be attributed to standard atomic 
HG modified by the fact that in clusters the atoms are 
disposed closer to each other while the physical origin 
of HG remains the same as in a gas jet. Only very re- 
cently, first experimental observations of the third har- 
monic (TH) generation from argon clusters subject to a 
strong laser field were reported in Ref. [TjJ where reso- 
nant enhancement of the TH yield, occurring when the 
Mie-frequency of the expanding cluster approaches three 
times the laser frequency, has been demonstrated using 
a pump-probe setup. The enhancement of the single- 
cluster response studied in theory before [H, 13|, 3 is, 
however, modified in the experiment by phase matching 
effects. The latter point complicates experimental stud- 
ies of nanoplasma radiation while in computations one 
can first examine the single-cluster response and may in- 
clude propagation effects in a second step. 

In this work, we concentrate on the single-cluster ra- 
diation in a short, intense laser pulse. The main ques- 
tion we address here is whether the nonlinear dynamics 
of laser-driven nanoplasmas lead to a substantial har- 
monic emission and which conditions are required to op- 
timize the signal. We study HG by three-dimensional 
particle-in-cell (PIC) simulations for large Arjv clusters 
(with the number of atoms N w 10 4 -10 5 and radii 
R Q w 6-10 nm) irradiated by linearly polarized, n = 8 
cycle sin 2 -laser pulses with an electric field E\(t) — 
E sin 2 (wii/2n) cos(a;i£) and wavelength A = 800 nm. 
Previous PIC simulations of HG from laser-driven clus- 
ters 0, [IH were performed for two-dimensional (rod-like) 
clusters. 

In a first step the laser field ionizes all neutral atoms 
to Ar + via over-the-barrier ionization once \E\(t)\ > 
Ip(Z)Z~ 1 /4 [lij], where Z is the charge number and 
I p (Z) is the respective ionization potential (atomic units 
are used). Since the laser field alters the charge equi- 
librium a space charge field E sc (r,t) is created. Ions of 
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FIG. 1: Normalized harmonic power P(ui) = \A(co)/N\ 2 
vs harmonic order for an Ari.7256 cluster at laser intensi- 
ties 2.5 x 10 14 Wcm -2 (HI), 2.5 x 10 15 Wcm~ 2 (H2), and 
7.5 x 10 17 Wcm~ 2 (H3). The spectrum H4 is for an Ar 920 96 
cluster at the intensity 7.5 x 10 17 Wcm~ 2 . 



higher charge states (Z > 1) are produced if the con- 
dition \E x (t) + E sc (R 3 ,t)\ > / 2 (Z)Z~7 4 is satisfied at 
the individual ion locations Rj. Since a PIC electron 
has the same charge to mass ratio as a real electron its 
equation of motion is fi = —E\(t) — E sc (ri,t). Each 
PIC ion of mass M and charge Zj(t) moves according 
MRj = Zj(t) [Ei(t) +E sc (Rj,t)]. The dipole approx- 
imation is applied, which is well justified for the laser 
wavelengths, cluster sizes and charge densities considered 
in this work. 

Because of azimuthal symmetry the total dipole ac- 
celeration A(t) is along the polarization direction of the 
laser pulse (x-axis). Its Fourier transformed amplitude 
A(u>) yields the dipole radiation power at the frequency 
u>. Figure [1] shows the normalized harmonic power 
P(ui) — \A{lo)/N\ 2 vs the harmonic order at various laser 
intensities. The spectra HI, H2, and H3 are the results 
for an Ari7256 cluster (i?o = 6.2 nm) while the spectrum 
H4 corresponds to an Arg 2 o96 cluster (R — 10.9 nm). 
At the intensity 2.5 x 10 14 Wcm~ 2 (HI) one observes a 
pronounced second harmonic. Upon increasing the laser 
intensity the third harmonic power also increases and be- 
comes comparable to the fundamental at higher laser 
intensities (H2). Increasing the intensity further, the 
fifth harmonic appears in the spectrum (H3). For the 
same intensity but a bigger cluster, the third harmonic 
is strongly suppressed as compared to the fifth harmonic 
(H4). 

The enhancements of particular harmonics are due to 
the resonance between their frequencies and the Mie- 
frequency of the expanding nanoplasma 12]. To prove 
this we retrieve the temporal information of the radiation 
by a time-frequency (TF) analysis of the spectra H1-H4 
in Fig. [H Figure [2K-d shows the respective TF diagrams. 
The scaled Mie- frequency u>yiie(t) /u>\ vs time is included 
in the plots. Note that a Mie-frequency can only be de- 




12345671234567 

Laser cycles Laser cycles 



FIG. 2: Time-frequency diagrams corresponding to the spec- 
tra HI - H4 of Fig[T]with the normalized power log 10 P. 

fined unambiguously as long as the ionic background re- 
mains homogeneously charged. From the simulations we 
find that the charge homogeneity is well satisfied within 
the initial cluster radius Rq while in the outer regions of 
the expanding cluster this is not the case. We therefore 
define WMic(i) = \/Qb{t) / Ro with Qb(t) the total ionic 
charge inside the sphere of radius Rq within which the 
cloud of well-bound electrons oscillates. 

At the intensity 2.5 x 10 14 Wcm~ 2 (Fig. [2^, correspond- 
ing to HI in Fig. [I]) the laser field yields only Ar + ions. 
This first ionization by the laser field alone gives rise to an 
abrupt jump of WMie(i)/wi up to the value tJMic/^i ~ 1-8, 
followed by a slower increase above this value due to ion- 
ization ignition [lit [l7| and finally ends with a plateau. 
The cluster expansion is so slow that tOMie{t) / to\ does not 
drop within the time interval plotted. The second har- 
monic power peaks when LOm e {t) / lo\ — 2 is met. This is 
because the electron dynamics contain oscillations both 
at the driving frequency and at the eigenfrcqucncy, de- 
pending on the initial conditions and the form of the 
potential. In long pulses such oscillations at the eigenfre- 
quency are damped due to various relaxation processes 
and do not contribute to the measured harmonic sig- 
nal. Figure [2Jd shows an enhanced third harmonic (H2 in 
Figd} at the time when ojyn e (t)/uj\ w 3. The third har- 
monic power starts increasing again around the 6th cycle 
when u>Mie(i) passes through the same resonance due to 
the cluster expansion. Figures^, d show the TF spectro- 
grams corresponding to spectra H3 and H4 (in FigLT]) at 
the intensity 7.5 x 10 17 Wcm~ 2 . Figured): clearly shows 
enhanced third and fifth harmonic emission (at 3-4 and 
6-7 cycles, respectively) when the scaled Mie-frequency 
approaches the respective odd numbers. For the bigger 
cluster Fig[2jl shows enhanced emission, preferentially 
following the Mie-frequency. As in the case with the sec- 
ond harmonic shown in Fig(2k- this is a consequence of 
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FIG. 3: TF spectrograms for two PIC electrons of the Ari7256 
cluster irradiated by the laser intensity 5 x 10 16 Wcm -2 . 
The dimensionless electron distance from the cluster center 
r(t)/Ro is indicated by the solid line. 



the undamped oscillations at the time-dependent eigen- 
frequency in the expanding cluster potential. However, 
pronounced emission starts around the 4th cycle when 
^Mie(i) meets 5u>i. Despite the same laser intensity as in 
Fig.[2j: u>Mic{t) does not reach 3u>i during the slower ex- 
pansion of the bigger cluster in Fig.[2ji. Note that the ful- 
fillment of the resonance condition u>Mie{t) — mw i (with 
m integer) is not sufficient for the emission of harmonics. 
A necessary condition is the presence of some nonlin- 
earity in the electron motion (i.e., anharmonicity in the 
effective potential). This nonlinearity may either origi- 
nate from electrons jutting out of the core during their 
motion, sensing the Coulomb tail [lij ]. or from the inho- 
mogeneous charge distribution within the ion core [l3T |. 

To clarify where the major contribution to the spec- 
trum comes from we consider the radiation of individual 
PIC particles. Figure [3] shows the TF analysis of the ra- 
diation from two PIC electrons. It is clearly seen that 
electrons radiate harmonics as long as they remain in- 
side the cluster. Leaving the cluster, they emit an in- 
tense flash with an almost continuous spectrum that ex- 
tends up to significantly higher frequencies than present 
in the net harmonic spectrum shown in FigQ] After liber- 
ation, the electrons emit only the fundamental frequency 
(linear Rayleigh scattering). A liberated electron may 
rescatter, giving rise to a second flash in the TF spectro- 
grams (FigO)). However, no indications of such intense 
flashes are visible in the net spectra shown in Figs. Q] and 
[2] The mechanism behind outer ionization (nonlinear 
resonance) allows to interpret these findings. Individ- 
ual electrons move in the selfconsistent field which could 
be subdivided into its slowly-varying and its oscillating 
component. The slowly- varying field is induced by the 
quasistatic part of the space charge. Ideally, it is a sta- 
tionary potential well if one ignores a slow evolution of 
the charge distribution due to the inner and outer ion- 
ization and the cluster expansion. The oscillating part 
of the field is a superposition of the incident laser field 
and the field induced by the oscillating electron cloud. 
In the low-frequency limit WMic(0 3> these two contri- 



butions are known to almost compensate each other so 
that the amplitude of the net oscillating field inside the 
cluster is small compared both to the applied laser field 
and the quasistatic space charge field. If the electron 
energy in the quasistatic well is far from the resonance 
its trajectory remains weakly disturbed by the oscillat- 
ing field. This causes HG with rapidly decreasing yield 
as a function of the harmonic order so that even the 7th 
and 9th harmonic are barely present in the spectra. As 
soon as the electron energy approaches the resonance, 
the same small perturbation results in a strong effect: the 
electron motion becomes strongly disturbed and stochas- 
tic. A stochastic near-resonance behavior is a well-known 
property of nonlinear systems driven by time-dependent 
forces [18]. In the cluster case it leads to almost prompt 
and irreversible outer ionization. Hence, passing through 
the resonance, an individual electron, upon leaving the 
cluster potential, emits radiation due to its strong accel- 
eration, seen as a flash in the TF spectrograms (Fig. [3]). 
However, exactly because of the stochastic nature of non- 
linear resonance the electrons' trajectories are very sen- 
sitive to the initial conditions with which the nonlinear 
resonance is entered. As a result, flashes from different 
electrons are incoherent (the corresponding amplitudes 
have nearly random phases), and, being added up in the 
total dipole acceleration, they disappear [3] . This shows 
that exactly the same mechanism behind efficient energy 
absorption by and outer ionization from clusters, namely 
nonlinear resonance 0, H, H, 0] > restricts HG from them 
by breaking the coherent electron motion once it becomes 
strongly anharmonic. Only well-bound electrons trapped 
inside the ionic core with energies far from the resonance 
contribute to the net, coherent radiation of the cluster. 

The time-dependent enhancements of particular har- 
monics analyzed above may be used for the reconstruc- 
tion of the maximum cluster charge density (i.e, max- 
imum Mie-frequency) . The TF spectrograms discussed 
above in Fig. [2] were obtained from our numerical ex- 
periments while it is not possible to record them in a 
real "single shot" experiment. To that end a more re- 
alistic pump-probe setup (see, e.g 
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20[ |) is envis- 
aged with the following PIC simulations. We combine 
the 800 nm near-infrared laser pulse (M-pulse, pump) 
with two UV pulses (X-pulses, probe) E V:Z (t — td) — 
Eo y ,z exp[— (ln2)(wi(t — td)/n) 2 ] sin[mwi(t — td)] of wave- 
length Xx = 800nm/m (with m = 3,5), applied at dif- 
ferent time delays td- To distinguish the contributions of 
M and X-pulses to the dipole acceleration we apply them 
polarized along y and z axis, respectively. The two X- 
pulses of different harmonic frequencies m ~ 3, 5 trigger 
independent Mie-oscillations of electrons in the respec- 
tive directions. The total power P VtZ (uj) = \A y , z (ui)/N\ 2 
is recorded at various time delays. Whenever Qb(t) is 
such that UMic{t) — muj\ for a certain time delay the cor- 
responding detector is expected to measure the maximum 
dipole radiation. Applying the two X-pulses simultane- 
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FIG. 4: Spectrogram of the dipole power log 10 P y 
y, z vs time delay td and frequency u/uj\ with the pump-pulse 
intensities (a,b) 2.5 x 10 15 Wcm~ 2 (cf. Fig©)), (c,d) 7.5 x 
10 17 Wcm~ 2 (cf. Fig[2};) and the probe-pulse intensity 2.5 x 
10 12 Wcm" 2 with m = 3 (a,c) and m = 5 (b,d). 



ously has the computational advantage of getting twice 
as much information in a single run than with just one 
probe pulse. Experimentally, a comb of low harmonics 
(e.g., from a gas jet) could be applied as a probe pulse. 

Figures 2k,b show log 10 P VtZ (uj) vs td and ui using 
X-pulses of frequencies 3uj\ and 5u>i, respectively. For 
lux = 3lji a resonance enhancement of the third har- 
monic power when u>Mi e /u)\ crosses ujx at id > 3.5 cycles 
is clearly seen, showing that the information revealed in 
the TF spectrogram of Fig. [^b is accessible in the more 
realistic pump-probe set-up. The enhancement is maxi- 
mum when cjMic/wi (i.e., the charge density p) is maxi- 
mum around td ~ 5 laser cycles. The frequency tux = 5u>i 
of the X-pulse along z in FigfJJ), being off-resonant with 
WMic, shows (an order of magnitude) less intense radia- 
tion. Figured^, allows the conclusion that the maximum 
cluster charge density was at least in the vicinity cor- 
responding to WMic ~ 3^1, i.e., max[p(t)] ~ 27ujf/4ir at 
that time. Figures5b,d are the analogues of Figs|Ui,b but 
at an M-pulse intensity 7.5 x 10 17 Wcm -2 , correspond- 
ing to Fig.[5fc. The enhanced radiation is clearly seen at 
those times when WMie/^i meets the X-pulse (along y) fre- 
quency Wx/ui = 3. However, the maximum of WMie/wi 
resonates with the other X-pulse (along z, FigHJi) fre- 
quency u>x = 5u>i at td ~ 3 cycles, and, indeed, the max- 
imum of the radiated power recorded in FigHJi is about 
an order of magnitude higher than in Fig.|4j:. 

In conclusion, we investigated numerically harmonic 
emission from laser-driven cluster nanoplasmas. The 
main contribution to the harmonic signal comes from 
electrons deeply bound inside the cluster potential. Such 
electrons emit only low-order harmonics with consider- 
able efficiency since they sense a weakly anharmonic 
potential only. In contrast, electrons passing through 



the nonlinear resonance and leaving the cluster move 
along strongly disturbed trajectories. Their radiation, 
although intense and broad in wavelength, is incoherent 
due to the stochastic nature of the nonlinear resonance 
and therefore does not contribute to the net signal from 
the whole cluster. The time-frequency analysis of the 
dipole acceleration shows enhanced third and fifth har- 
monic emission when the Mie-frequency of the expand- 
ing cluster meets the respective harmonic frequencies, 
which is consistent with both previous theoretical studies 
llj 13 j id | and recent experiment 11 1. A pump-probe 
experiment is proposed to measure the cluster charge 
density by detecting the dipole radiation at different time 
delays. This method can be used to monitor the in- 
ner and outer ionization dynamics of clusters. A control 
scheme to achieve the highest possible charge states by 
driving the cluster as resonantly as possible is currently 
investigated. 
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